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Abstract This work deals with a computational in- 
vestigation on the small ballistic reentry Brazilian ve- 
hicle SARA (acronyms for SAtelite de Reentrada At- 
mosfcrica). Hypersonic flows over the vehicle SARA at 
zero-degree angle of attack in a chemical equilibrium 
and thermal non-equilibrium are modeled by the Di- 
rect Simulation Monte Carlo (DSMC) method, which 
has become the main technique for studying complex 
multidimensional rarefied flows, and that properly ac- 
counts for the non-equilibrium aspects of the flows. The 
emphasis of this paper is to examine the behavior of the 
primary properties during the high altitude portion of 
SARA reentry. In this way, velocity, density, pressure 
and temperature field are investigated for altitudes of 
100, 95, 90, 85 and 80 km. In addition, comparisons 
based on geometry are made between axisymmetric and 
planar two-dimensional configurations. Some significant 
differences between these configurations were noted on 
the flowfield structure in the reentry trajectory. The 
analysis showed that the flow disturbances have differ- 
ent influence on velocity, density, pressure and temper- 
ature along the stagnation streamline ahead of the cap- 
sule nose. It was found that the stagnation region is a 
thermally stressed zone. It was also found that the stag- 
nation region is a zone of strong compression, high wall 
pressure. Wall pressure distributions are compared with 
those of available experimental data and good agree- 
ment is found along the spherical nose for the altitude 
range investigated. 
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1 Introduction 
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The microgravity field has become nowadays one of 
the most promising new areas for the commercializa- 
tion of space. The development of scientific and tech- 
nological experiments carried out under conditions of 
reduced gravity has been proposed by the Brazilian 
Space Agency by employing a recoverable orbital sys- 
tem. The orbital system, built in a platform with a cap- 
sule shape, is designed in order to stay in orbit during 
the time needed for the execution of the experiments. 
After that, the capsule is sent back to the Earth and 
recovered. 

The aerothermodynamic aspects during the ballistic 
reentry flight, in the development of a capsule platform, 
offer exciting challenges to the aerodynamicists. In the 
Earth atmosphere reentry, the capsule undergoes not 
only different velocity regimes - hypersonic, supersonic 
and subsonic - but also different flow regimes - free 
molecular flow, transition and continuum - and flight 
conditions that may difficult their aerodynamic design. 
Therefore, the capsule aerodynamic design, of great im- 
portance for the flight through several atmosphere lay- 
ers, has to consider important aspects such as stabil- 
ity, drag and heating loads. In this context, a combi- 
nation of engineering tools, experimental analysis, and 
numerical methods is used in the design of high altitude 
reentry capsule aerodynamics. Nevertheless, due to dif- 
ficulties and high costs associated to the experimental 
work at high speed flows, a numerical analysis becomes 
imperative. 
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Many experimental and numerical studies [1,2,3,4, 
5,6,7,8,9, 10, 11, 12] have been dedicated to the aerother- 
modynamic of vehicles reentering the Earth atmosphere. 
Nevertheless, for the particular case of SARA capsule, 
only a few studies are available in the current litera- 
ture [13,14,15,16,17,18]. For the purpose of this intro- 
duction, it will be sufficient to describe only some of 
these studies. 

Sharipov [13] has investigated the flowfield struc- 
ture over the SARA capsule by employing the DSMC 
method. It was assumed the monatomic gas argon as 
the working fluid and freestream Mach numbers of 5, 
10, and 20. Even considering that the real gas effect in a 
reentry capsule is not well represented by a monotonic 
gas, that investigation might be considered as the first 
contribution to the aerothcrmodynamic analysis of the 
SARA capsule in high altitudes. 

Pimcntel et al. [14] have performed inviscid hyper- 
sonic flow simulations over the SARA capsule modeled 
by employing the planar two-dimensional (2-D) and the 
axisymmctric Eulcr equations. Results were presented 
for an altitude of 80 km, Mach numbers of 15 and 18, 
and angle of attack of and 10 degrees. They also con- 
sidered air as working fluid composed of five species 
(N2, O2, O, N, and NO) with their reactions of disso- 
ciation and recombination. Pressure and temperature 
contour maps were presented for 2-D and axisymmet- 
ric flow. 

Finally, by using axisymmetric Navier-Stokes equa- 
tions, Tchuen et al. [15] have investigated the flowfield 
structure over the SARA capsule by considering hy- 
personic flow at zero angle of attack in chemical and 
thermal non-equilibrium. It was assumed air as work- 
ing fluid composed of seven species (N 2 , O2, O, N, NO, 
NO + , and e~) associated with their reactions of dis- 
sociation and recombination. Results for pressure, skin 
friction, and heat transfer coefficients were presented 
for different combinations of Mach numbers of 10, 20 
and 25 with altitudes of 75 and 80 km. 

The interest in the majority of these SARA studies 
has basically gone into considering the flow in the con- 
tinuum flow regime. Conversely, in the transition flow 
regime, there is a little understanding of the physical 
aspects of a hypersonic flow past to the SARA capsule 
related to the severe aerothcrmodynamic environment 
in the reentry trajectory. In this fashion, the purpose of 
the present account is to investigate the flowfield struc- 
ture of a hypersonic flow over the SARA capsule in the 
transition flow regime, i.e., between the free collision 
flow and the continuum flow regime. The primary goal 
is to assess the sensitivity of the primitive properties, 
such as, velocity, density, pressure, and temperature 
due to changes on the altitude representative of a typ- 



ically reentry trajectory of the SARA capsule. There- 
fore, the present study focuses on the low-density region 
in the upper atmosphere, where the non-equilibrium 
conditions are such that the traditional Computational 
Fluid Dynamics (CFD) calculations arc inappropriate 
to yield accurate results. In such a circumstance, the 
Direct Simulation Monte Carlo (DSMC) method will 
be employed to calculate the planar two-dimensional 
and the axisymmetric hypersonic flow over the SARA 
capsule. 

It is important to mention that, for the particular 
case of SARA, this is the first work to present the flow- 
field structure in the transition flow regime by consid- 
ering two species (N2 and O2) as working fluid, and 
rotation and vibration internal modes of energy. 

2 Geometry Definition 

The SARA reentry capsule is an axisymmetric design 
consisting of a spherical nose with a 11.4-degree half- 
angle conical afterbody. The nose radius R is 0.2678 




h 

L 



Fig. 1 Drawing illustrating a schematic view of the capsule 
(top), and the important parameters (bottom). 
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m, the afterbody base has a radius Rb of 0.5035 m, 
and the total length L is 1.410 m. Figure 1 illustrates 
schematically the capsule shape and the main impor- 
tant physical and geometric parameters related to the 
hypersonic flow on the capsule. The main physical pa- 
rameters are defined as follows: is the freestream 
Mach number, Kn^ stands for the Knudsen number, 
a is the angle of attack, T w is the wall temperature, 
and finally, a n and at represent the parameters related 
to the gas-surface interaction. In this fashion, the flow- 
field structure around the capsule may be affected by 
the effects of compressibility, rarefaction, gas-surface in- 
teraction, etc. 

3 Freestream and Flow Conditions 

Flow conditions represent those experienced by the cap- 
sule in the reentry trajectory from 100 to 80 km of alti- 
tude. This range of altitude is associated with the tran- 
sition flow regime, which is characterized by the overall 
Knudsen number the order of or larger than 10 -2 . 

Freestream flow conditions used for the numerical 
simulation of flow past the capsule are those given by 
Hirschel [19] and summarized in Tables 1 and 2, and 
the gas properties [20] are shown in Table 3. Referring 
to this set of tables, Too, Poo, Poo, noo, Aqo, and 
stand respectively for temperature, pressure, density, 
number density, molecular mean free path, and veloc- 
ity, and x, m , d an d w account respectively for mass 
fraction, molecular mass, molecular diameter and vis- 
cosity index. 

The velocity-altitude map for the SARA capsule [17, 
21] is demonstrated in Fig. 2. This velocity-altitude 
map was generated based on predefined conditions, such 

Table 1 Freestream flow conditions [19]. 

Altitude (km) T X {K) Poo (N/m 2 ) Poo (kg/m 3 ) 



80 


180.7 


1.03659 


1.999 x 10- 


5 


85 


180.7 


0.41249 


7.956 x 10- 


6 


90 


180.7 


0.16438 


3.171 x 10- 


6 


95 


195.5 


0.06801 


1.212 x 10- 


6 


100 


210.0 


0.03007 


4.989 x 10- 
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Table 2 Freestream flow conditions [19] (cont'd). 



Altitude (km) 


iioo(m 3 ) 


Aoo(m) 




Voo(m/s) 


80 


4.1562 x 10 21 


3.085 x 10- 


3 


7820 


85 


1.6539 x 10 20 


7.751 x 10- 


3 


7864 


90 


6.5908 x 10 19 


1.945 x 10- 


2 


7864 


95 


2.5197 x 10 19 


5.088 x 10- 


2 


7866 


100 


1.0372 x 10 19 


1.236 x 10- 


1 


7862 



Table 3 Gas properties [20] 





X 


m (kg) 


d (m) 


tc 


o 2 


0.237 


5.312 x 10" 26 


4.01 x 10" 10 


0.77 


N 2 


0.763 


4.650 x 10" 26 


4.11 x 10" 10 


0.74 



as an initial velocity of 7626.30 m/s, an initial altitude 
of 300 km, a drag coefficient of 0.80, a lift coefficient of 
0.0, a reference area of 0.785 m 2 , a total mass of 150 
kg, and finally, a gravity acceleration of 9.81 m/s 2 . 

According to Fig. 2, for altitudes of 100, 95, 90, 85, 
and 80 km, the freestream velocity Vx is 7862, 7866, 
7864, 7864, and 7820 m/s, respectively. These values 
correspond to a freestream Mach number of 27.1, 
28.1, 29.2, 29.2, and 29.0, respectively. In the present 
account, the capsule surface was kept at a constant wall 
temperature T w of 800 K for all cases investigated. This 
temperature is chosen to be representative of the sur- 
face temperature near the stagnation point of a reentry 
capsule. According to Machado and Boas [22], in the 
stagnation region of the SARA capsule, temperature 
may reach a maximum value around 950 K. 

The overall Knudsen number Kn is defined as the 
ratio of the molecular mean free path A in the freestream 
gas to a characteristic dimension of the flowficld. In the 
present study, the characteristic dimension was defined 
as being the nose radius R. For the altitudes investi- 
gated, 100, 95, 90, 85, and 80 km, the overall Knudsen 
numbers correspond to Kn R of 0.4615, 0.1899, 0.0726, 
0.0289, and 0.0115, respectively. Finally, the Reynolds 
number Re R correspond to 92, 224, 609, 3442 and 15249 
for altitudes of 100, 95, 90, 85 and 80 km, respectively, 
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Fig. 2 The velocity-altitude map for the SARA capsule. 
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based on conditions in the undisturbed stream with the 
nose radius R as the characteristic length. 



4 Computational Method and Procedure 

The governing equation in the transition flow regime 
is the Boltzmann equation [23]. The Boltzmann equa- 
tion is a nonlinear integral-differential equation, with 
one dependent variable given by the particle distribu- 
tion function. The particle properties are determined in 
a six-dimensional space, named the phase space, com- 
posed by three dimensions for particle coordinates, and 
three dimensions for particle velocities. 

In order to circumvent the difficulty of a direct solu- 
tion of the Boltzmann equation, the Direct Simulation 
Monte Carlo (DSMC) method is one of the alternative 
approaches for solving the Boltzmann equation by sim- 
ulating the behavior of individual particles. The DSMC 
method, pioneered by Bird [20], has become the appro- 
priate choice for problems involving complex multidi- 
mensional flows of rarefied hypersonic aerothermody- 
namics. In addition, CFD methods, which rely on con- 
tinuum relations to compute the flowfield structure, will 
not provide accurate results for these flows in the upper 
atmosphere, since the assumptions made in developing 
the differential equations, on which CFD methods are 
based, break down on rarefied conditions. 

In the DSMC method, the gas is modeled at the 
microscopic level by simulated particles. Each simu- 
lated particle represents a very large number of real 
molecules or atoms. These representative molecules are 
tracked as they move, collide and undergo boundary 
interactions in simulated physical space. The molecular 
motion, which is considered to be deterministic, and 
the intermolecular collisions, which are considered to 
be stochastic, are uncoupled over the small time step 
used to advance the simulation and computed sequen- 
tially. 

In the present account, collisions are modeled by us- 
ing the variable hard sphere (VHS) molecular model [24] 
and the no time counter (NTC) collision sampling tech- 
nique [25]. The VHS model employs the simple hard 
sphere angular scattering law. Therefore, all directions 
are equally possible for post-collision velocity in the 
center-of-mass frame of reference. Nevertheless, the col- 
lision cross section depends on the relative speed of col- 
liding molecules according to some power law. The ex- 
ponent is calculated by matching the viscosity of the 
simulated gas to that of its real counterpart. 

The mechanics of the energy exchange processes be- 
tween kinetic and internal modes for rotation and vibra- 
tion are controlled by the Borgnakke-Larsen statistical 



model [26]. The essential characteristic of this model 
is that a fraction of collisions is treated as completely 
inelastic, and the remainder of the molecular collisions 
is regarded as elastic. For a given collision, the proba- 
bilities are designated by the inverse of the relaxation 
numbers, which correspond to the number of collisions 
necessary, on average, for a molecule to relax. In this 
study, the relaxation numbers of 5 and 50 were used for 
the rotation and vibration, respectively. Simulations are 
performed using a non-reacting gas model consisting of 
only two chemical species, N 2 and O2. 

Finally, the frccstrcam coefficient of viscosity fi^, 
and the freestream mean free path Aoo are evaluated 
from a consistent definition [20,27] by using the VHS 
collision model with temperature exponents (viscosity 
index) of 0.74 and 0.77 for N 2 and O2, respectively. 



5 Computational Flow Domain and Grid 

For the numerical treatment of the problem, the flow- 
field around the capsule is divided into an arbitrary 
number of regions, which are subdivided into computa- 
tional cells. The cells are further subdivided into sub- 
cells, two subcells/cell in each coordinate direction. The 
cell provides a convenient reference for the sampling of 
the macroscopic gas properties. On the other hand, the 
collision partners are selected from the same subcell 
in order to establish the collision rate. Therefore, the 
physical space network is used to facilitate the choice 
of molecules for collisions as well as for the sampling of 
the macroscopic flow properties, such as temperature, 
pressure, density, etc. 

The computational domain used for the calculation 
is made large enough so that the capsule disturbances 
do not reach the upstream and side boundaries, where 
freestream conditions are specified. In this manner, the 
computational domain changed according to the rar- 
efaction degree of the flow on the capsule. A schematic 
view of the computational domain is depicted in Fig. 3. 
Advantage of the flow symmetry is taken into account, 
and molecular simulation is applied to one-half of a full 
configuration. According to this plot, side I is defined 
by the capsule surface. Diffuse reflection with complete 
thermal accommodation is the condition applied to this 
side. In a diffuse reflection model, the molecules are re- 
flected equally in all directions. In addition, the final 
velocity of the molecules is randomly assigned accord- 
ing to a half-range Maxwellian distribution based on 
the wall temperature. Side II is a plane of symmetry, 
where all flow gradients normal to the plane are zero. At 
the molecular level, this plane is equivalent to a spec- 
ular reflecting boundary. Side III is the freestream side 
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Fig. 3 Drawing illustrating the computational domain. 

through which simulated molecules enter and exit. Fi- 
nally, the flow at the downstream outflow boundary, 
side IV, is predominantly supersonic and vacuum con- 
dition was assumed at this boundary [28]. As a result, 
simulated molecules can only exit at this boundary. 

The numerical accuracy in DSMC method depends 
on the cell size chosen, on the time step as well as 
on the number of particles per computational cell. In 
the DSMC algorithm, the linear dimensions of the cells 
should be small in comparison with the length scale of 
the macroscopic flow gradients normal to streamwise 
directions. Therefore, the cell dimensions should be the 
order of or smaller than the local mean free path [29, 
30]. Furthermore, the time step should be chosen to 
be sufficiently small in comparison with the local mean 
collision time [31,32]. In general, the total simulation 
time, discretized into time steps, is based on the phys- 
ical time of the real flow. Finally, the number of simu- 
lated particles has to be large enough to make statistical 
correlations between particles insignificant. 

These effects were investigated in order to determine 
the number of cells and the number of particles required 
to achieve grid independent solutions. The grid gener- 
ation scheme used in this study follows that procedure 
presented by Bird [20]. Along the body surface (side I) 
and the outer boundary (side III), point distributions 
are generated in such way that the number of points on 
each side is the same; ^-direction in Fig. 3. Then, the 
cell structure is defined by joining the corresponding 
points on each side by straight lines and then divid- 
ing each of these lines into segments which are joined 
to form the system of quadrilateral cells; 77-direction in 
Fig. 3. The distribution can be controlled by a number 
of different distribution functions that allow the concen- 



tration of points in regions where high flow gradients or 
small mean free paths are expected. Grid independence 
was tested by running the calculations with a coarse, 
standard and a fine grid. A discussion of grid effects on 
the aerodynamic surface quantities is described in the 
Appendix. 

6 Computational Results and Discussion 

This section focuses on the effects that take place in the 
primary properties due to variations on the altitude of 
the capsule. Primary properties of particular interest in 
the transition flow regime are velocity, density, pressure 
and temperature. In this scenario, this section discusses 
and compares differences in these quantities due to rar- 
efaction effects. 

6.1 Velocity Field 

Normal velocity profiles along the stagnation stream- 
line and their dependence on rarefaction are illustrated 
in Figs. 4(a-c) for altitudes of 100, 90 and 80 km, re- 
spectively. In this set of plots, the normal velocity v is 
normalized by the freestream velocity Voc , and the dis- 
tance x along the stagnation streamline is normalized 
by the nose radius R. In addition, for comparison pur- 
pose, the velocity profiles for the 2-D geometry are also 
included in the plots. Also, in order to emphasize points 
of interest, a different scale is used in the abscissa axis. 
It is important to note that Voa is slightly different for 
each altitude (see Tab. 1) and, therefore, the compari- 
son is made in terms of ratio. Also, the velocity profiles 
for the 95 km and 85 km cases are intermediate to the 
other cases and, therefore, they will not be shown. 

According to these plots, it is seen that the rarefac- 
tion effect influences the flowfield far upstream. The ex- 
tent of this effect decreases with decreasing the Knud- 
sen number Kna, i.e., by decreasing the altitude. As the 
altitude decreases, particles reflecting from the capsule 
surface diffuse less upstream due to the high density 
at the vicinity of the capsule nose, as will be seen sub- 
sequently. As a result, the extent of the flowfield dis- 
turbance due to the presence of the capsule decreases 
and becomes smaller in terms of the body dimension R. 
It should be mentioned in this context that, far from 
the capsule nose, the velocity ratio v/Voo is equal to 
one, and at the stagnation point equal to zero. The 
region defined from the section in which the velocity 
ratio starts decreasing from 1, say v/Vao = 0.99, up 
to the stagnation point, v/Vrx, = 0, is defined herein 
as the upstream disturbance region, i.e., the region af- 
fected by the presence of the capsule. For instance, for 
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Fig. 4 Normal velocity (— v/Voo) profiles along the stagna- 
tion streamline for altitudes of (a) 100 km, (b) 90 km, and 
(c) 80 km. 



Fig. 5 Tangential velocity (u/Voo) profiles along the capsule 
surface for sections corresponding to s/R of (a) 0.4, (b) 0.8 
and (c) 1.6. 



the axisymmetric geometry, the extent of the upstream 
disturbance, based in a velocity reduction of 1% (v/Voo 
= 0.99), corresponds to sections x/R of 1.86, 1.20, 0.75, 
0.45, and 0.26 for altitudes of 100, 95, 90, 85, and 80 
km, respectively. Nevertheless, for the 2-D geometry, it 
changes to around 3.51, 2.12, 1.21, 0.71, and 0.45 for al- 
titudes of 100, 95, 90, 85, and 80 km, respectively. Con- 



sequently, the extension of the upstream disturbance is 
more pronounced for a 2-D flow than for an axisym- 
metric flow. It should be remarked that this effect is 
less pronounced for an axisymmetric geometry because 
of the "relieving effect" , which is a characteristic of all 
three-dimensional flows. For the flow over an axisym- 
metric capsule, the "addition" of a third dimension pro- 
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vidcs the molecules with extra space to move through it. 
In contrast, for the 2-D geometry, the molecules move 
only around the top or the bottom surface. 

The outer extent of the flowheld disturbance over 
the capsule surface is demonstrated in Figs. 5(a-c) as 
a function of the altitude. In this set of plots, the tan- 
gential velocity u is normalized by the freestream ve- 
locity Voo, and the height rj in the off-body direction 
(77-direction in Fig. 3(a)) is normalized by the nose ra- 
dius R. In an effort to emphasize points of interest, this 
set of plots presents data for three stations along the 
body surface defined by s/R of 0.4, 0.8, and 1.6, where s 
is the arc length along the body surface measured from 
the stagnation point. The first two sections are located 
on the spherical nose, and the last one on the afterbody 
surface. 

Interesting features can be drawn from this set of 
tangential velocity profiles. As the body slope decreases, 
the tangential velocity adjacent to the body surface in- 
creases. This is to be expected since the flow experi- 
ences an expansion as it moves downstream along the 
capsule surface. It is observed that the flow accelerates 
faster along the surface with the altitude rise. It should 
be noted that the tangential velocity Moo, defined as 
7] — > 00, is represented by the dashed line and shown 
for each station. Because of the body curvature, Uoo ( = 
Voo cos 9) varies as a function of the body slope. 

Another interesting characteristic in these plots is 
the similarity of the velocity profiles along the body 
surface. This is an indication that the velocity profiles 
may be expressed in terms of functions that, in appro- 
priate coordinates, may be independent of one of the 
coordinate directions. However, no attempts have been 
done to find such functions. 



6.2 Density Field 

The impact of rarefaction on density profiles along the 
stagnation streamline is displayed in Figs. 6(a-c) for 
altitudes of 100, 90 and 80 km, respectively. In this set 
of figures, density p is normalized by the freestream 
density poo. Again, as a base of comparison, density 
profiles for the 2-D geometry are also included in the 
figures. 

The predictions of density for the altitude range in- 
vestigated basically show no sign of a discrete shock 
wave. Except for the 80 km case, it is clearly seen a 
continuous rise in density from the freestream to the 
nose of the capsule, an indication of the diffuse nature 
of the shock wave, which corresponds to a character- 
istic of highly rarefied flows. By decreasing the alti- 
tude, both the freestream gas density and temperature 



change. As a result, the molecules interact much more 
with each other and collisions among them are more 
frequent. Therefore, density increases and the profile 
becomes steeper near the stagnation point. In addition, 
with the density rise near the stagnation point, the lo- 
cal mean free path decreases resulting in a lower local 
Knudscn number. Moreover, it is noticed that, unlike 
normal velocity, density has little effect on the extent 
of the upstream disturbance caused by the presence of 
the capsule, for the cases investigated. Much of the 
density increase in the shock layer occurs in a region 
smaller than one nose radius R, and after the temper- 
ature has reached its postshock value, as will be seen 
subsequently. 

Still referring to Figs. 6(a-c), it can be recognized 
that density rises to well above the continuum invis- 
cid limit for the cases investigated. As a point of refer- 
ence, the Rankine-Hugoniot relations give a postshock 
density that corresponds to the ratio p/poo of 5.96 for 
freestream Mach number of 27. Near the stagnation 
point, x/R ~ 0, a substantial density increase occurs, 
which is a characteristic of a cold- wall entry flow. In a 
typical entry flow, the body surface temperature is low 
compared to the stagnation temperature. This leads to 
a steep density gradient near to the body surface. For 
the present simulation, the ratio of wall temperature to 
stagnation temperature is around 0.033, which corre- 
sponds to a cold- wall flow. 

The deceleration of the freestream air molecules leads 
to an increase in the internal energy, such as trans- 
lational, rotational and vibrational modes. Since the 
vibrational phenomenon absorbs energy, the temper- 
ature does not increase to its perfect gas value, and 
consequently, the density rise across the shock wave is 
greater than that for a perfect gas flow at the same 
freestream Mach number. For air behaving as a per- 
fect gas at hypersonic speeds p/poc — > 6, while for 
atmospheric hypersonic flight p/ Poo is the order of 10 
to 14 or more [33] , because of the elevated temperatures 
and resultant activation of more degrees of freedom. As 
a reference, Cuda Jr. and Moss [34] have reported a 
simulation on hypersonic flow over blunt wedges at an 
altitude of 70 km, and they found a value for the stag- 
nation density two orders of magnitude greater than 
the freestream value. 

Variation of local density along the body normal di- 
rection, expressed as a ratio to the freestream value, 
is depicted in Figs. 7(a-c) as a function of the altitude. 
This group of diagrams presents data at three afterbody 
stations that correspond to the dimensionless arc length 
s/R of 0.4, 0.8, and 1.6. According to these diagrams, it 
is clearly noticed that the density also experiences sig- 
nificant changes in the direction perpendicular to the 
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Fig. 6 Density ratio (p/poo) profiles along the stagnation 
streamline for altitudes of (a) f 00 km, (b) 90 km, and (c) 80 
km. 



Fig. 7 Density ratio (p/poa) profiles along the capsule sur- 
face for sections corresponding to s/R of (a) 0.4, (b) 0.8 and 
(c) 1.6. 



wall as the flow moves downstream along the capsule 
surface. In the direction perpendicular to the wall, and 
for stations close to the stagnation region, the density 
is high adjacent to the wall and rapidly decreases in- 
side a layer of thickness smaller than one nose radius 
R, where the density approaches the freestream value 
for the altitudes investigated. This characteristic is ob- 



served when the body surface is very much colder than 
the stagnation temperature of the oncoming gas. As a 
result, the gas near the body tends to be much denser 
and cooler than the gas in the rest of the boundary 
layer. Conversely, for a station far from the stagnation 
region, s/R = 1.6, the density ratio is maximum at the 
wall, it decreases in the off-body direction r\ and then it 
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increases to a peak value inside the shock wave. After- 
wards, the density decreases dramatically and reaches 
the freestream density value as r\/R — > oo. 

Still referring to Fig. 7, it is very encouraging to 
observe that density is affected by changes in the al- 
titude, as would be expected. According to Fig. 7(a), 
which corresponds to the section s/R = 0.4, the density 
variation is in excess of one order of magnitude as com- 
pared to the freestream density for the altitude inves- 
tigated. In this region, close to the stagnation region, 
the compression combined with a relatively cool wall 
produces a maximum density ratio p/poo that lies in 
the range from 40 to 90. Because of the flow expansion 
along the afterbody surface, the density ratio adjacent 
to the surface decreases dramatically to the range of 6 
to 8 for the section corresponding to s/R = 1.6. 

6.3 Pressure Field 

Part of the large amount of kinetic energy present in 
a hypersonic freestream is converted by molecular col- 
lisions into high thermal energy surrounding the body 
and by flow work into increased pressure. In this man- 
ner, the stagnation line is a zone of strong compres- 
sion, where pressure increases from the freestream to 
the stagnation point due to the shock wave that forms 
ahead of the capsule. 

Pressure profiles along the stagnation streamline are 
depicted in Figs. 8(a-c) as a function of the altitude. 
In this set of plots, pressure p is normalized by the 
freestream pressure p^ . It may be recognized from these 
plots that basically there is a continuous rise in pressure 
from the freestream up to the stagnation point where 
the maximum value is attained. Near the stagnation 
point, a substantial pressure increase occurs with de- 
creasing the altitude, i.e., with decreasing the Knudsen 
number Khr. It is apparent from these plots that the 
general shape of the pressure distribution profiles is pre- 
served when the altitude decreases from 100 km to 80 
km. 

The extent of the upstream flowficld disturbance for 
pressure is significantly different from that presented by 
velocity and density. The domain of influence for pres- 
sure is higher than that for velocity and density, and 
lower than that presented for temperature. Similar to 
the density, much of the pressure increase in the shock 
layer occurs after the translational kinetic temperature 
has reached its postshock value, as will be shown sub- 
sequently. 

Local pressure, expressed as a ratio to the freestream 
value, for three stations located on the capsule surface 
is demonstrated in Figs. 9(a-c) as a function of the al- 
titude. It is apparent from these profiles that pressure 



1200 



1000 



(2 600 













































ft 




100 km 




-• — Axi. 
■ 2-D 

1 







1200 



.P. 800 



P 400 



200 




10 



Dimensionless Length (-x/R) 




10° 10 1 10 2 10 3 

Dimensionless Length (-x/R) 



lO 4 




Dimensionless Length (-x/R) 

Fig. 8 Pressure ratio (p/poa) profiles along the stagnation 
streamline for altitudes of (a) 100 km, (b) 90 km, and (c) 80 
km. 



is affected with decreasing altitude, as was mentioned 
earlier. For the station corresponding to s/R = 0.4, 
Fig. 9(a), the pressure variation is almost three orders 
of magnitude larger than the freestream pressure for 
an altitude of 80 km. In this region, at the vicinity of 
the stagnation point, the compression produces a max- 
imum pressure that is around 900 times the freestream 
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Fig. 9 Pressure ratio (p/poo) profiles along the capsule sur- 
face for sections corresponding to s/fl of (a) 0.4, (b) 0.8 and 
(c) 1.6. 



pressure for an altitude of 80 km. Due to the flow ex- 
pansion along the body surface, the pressure adjacent 
to the surface decreases significantly. For the particu- 
lar case of 80 km, the pressure ratio p/p x decreases to 
around 90 for the station corresponding to s/R = 1.6, 
as shown in Fig. 9(c), a reduction of basically one order 
of magnitude from station s/R of 0.4 to 1.6. 



6.4 Temperature Field 

The strong shock wave that forms ahead of the capsule 
at hypersonic flow converts part of the kinetic energy of 
the freestream air molecules into thermal energy. This 
thermal energy downstream of the shock wave is par- 
titioned into increasing the translational kinetic energy 
of the air molecules, and into exciting other molecular 
energy states such as rotation and vibration. 

Representative kinetic temperature profiles along the 
stagnation streamline are demonstrated in Figs. 10 and 
11 for 100, 90 and 80 km of altitude. Figures 10(a- 
c) exhibit the temperature distribution with a linear 
scale in the abscissa in order to emphasize the exten- 
sion of the upstream disturbance, while Figs, ll(a-c) 
depict the same temperature distribution in a loga- 
rithm scale in order to emphasize the temperature be- 
havior at the vicinity of the stagnation point. In addi- 
tion, in this set of pictures, temperature ratio stands for 
the translational temperature Tx, rotational tempera- 
ture Tr, vibrational temperature Ty or overall tem- 
perature To normalized by the freestream temperature 
r^. Also, filled and empty symbols correspond to tem- 
perature distributions for axisymmctric and 2-D geome- 
tries, respectively. It is apparent from these figures that 
thermodynamic non-equilibrium occurs throughout the 
shock layer, as shown by the lack of equilibrium of the 
translational and internal kinetic temperatures. Ther- 
mal non-equilibrium occurs when the temperatures as- 
sociated with the translational, rotational, and vibra- 
tional modes of a polyatomic gas are different from each 
other. In this context, an overall kinetic temperature To 
is defined [20] for a non-equilibrium g RS MS the weighted 
mean of the translational and internal temperature by 
the following expression, 

CtTt + (rTr + CvTy 
° = ? X7 Z7 W 

were ( is the degree of freedom and subscript T, R and 
V stand for translation, rotation and vibration, respec- 
tively. 

The overall kinetic temperature To is equivalent to 
the thermodynamic temperature only under thermal 
equilibrium conditions. It should be emphasized that 
the ideal gas equation of state does not apply to this 
temperature in a non-equilibrium situation. 

Referring to Figs. 10 and 11, in the undisturbed 
freestream far from the capsule, the translational and 
internal temperatures have the same value and are equal 
to the thermodynamic temperature. Approaching the 
nose of the capsule, the translational temperature rises 
to well above the rotational and vibrational tempera- 
tures and reaches a maximum value that is a function of 
the altitude. Since a large number of collisions is needed 
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Fig. 10 kinetic temperature ratio (T/Too) profiles along the 
stagnation streamline in a linear scale for altitudes of (a) 100 
km, (b) 90 km, and (c) 80 km. 



Fig. 11 kinetic temperature ratio (T/Too) profiles along the 
stagnation streamline in a logarithm scale for altitudes of (a) 
100 km, (b) 90 km, and (c) 80 km. 



to excite the molecules in the vibrational mode, from 
the ground state to the upper state, then the vibra- 
tional temperature increases much more slowly than ro- 
tational temperature. Still further downstream toward 
the nose of the capsule, the translational temperature 
decreases and reaches a value on the wall that is above 
the wall temperature (T^/Tqo « 4) for the 100 km and 



90 km of altitude, resulting in a temperature jump. 
Therefore, for x/R ~ 0, one has Tt ^ T w . In con- 
trast, for the 80 km case, the translation, rotation and 
vibrational temperature are basically the same at the 
vicinity of the nose, indicating the the thermodynamic 
equilibrium is achieved. 
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The substantial rise in the kinetic translational tem- 
perature occurred before the density rise, as shown in 
Figs. 6(a-c). For instance, the translational temperature 
ratio, Tt/Tqo, for the axisymmetric case reaches the 
maximum value of 119.3, 169.2, and 186.0, for altitude 
of 100, 90 and 80 km, respectively, around a distance 
of 0.40-R, 0.18i?, and 0.12i?, respectively, from the nose 
of the capsule, while the density ratio p/poo, for the 
same stations, is around 2.1, 2.3, and 2.6 for altitude 
of 100, 90 and 80 km, respectively. The translational 
kinetic temperature rise for a blunt body results from 
the essentially bimodal velocity distribution: the molec- 
ular sample consisting of mostly undisturbed freestream 
molecules with the molecules that have been affected by 
the shock and reflected from the capsule nose. In this 
scenario, the translational kinetic temperature rise is 
a consequence of the large velocity separation between 
these two classes of molecules. 

Particular attention is paid to the kinetic transla- 
tional temperature in the shock layer. In this respect, 
the kinetic translational temperature variation is taken 
normal to the capsule surface at stations corresponding 
to the dimcnsionlcss arc length s/R of 0.4, 0.8, and 1.6. 
Figures 12(a-c) depict the kinetic translational temper- 
ature profiles at the considered positions normal to the 
capsule surface along the 77-axis for the altitudes inves- 
tigated. According to these figures, it is observed that 
the downstream evolution of the flow displays a smear- 
ing tendency of the shock wave due to the displace- 
ment of the maximum value for the kinetic translational 
temperature. Also, it may be recognized from the tem- 
perature distribution in Figs. 12(a-c) that significant 
changes in the translational temperature profiles along 
the spherical nose occur within a layer adjacent to the 
capsule surface of a few nose radius R for the altitude 
range investigated. 

In order to bring out important features of the rar- 
efaction effects, particular attention is paid to the over- 
all kinetic temperature at the vicinity of the capsule 
nose. In this scenario, overall kinetic temperature con- 
tours, normalized by the freestream temperature Too, 
are plotted in Figs. 13(a-c) for altitude of 100, 90 and 
80 km, respectively. In this set of plots, the length x 
and height y are normalized by the nose radius R. In 
addition, the upper half part of the plots represents the 
axisymmetric case and the lower half part stands for the 
2-D case. Also, in order to emphasize points of interest, 
a different scale is used in the abscissa and ordinate axis 
of these plots. This set of plots clearly illustrates the dif- 
ferences in the temperature distribution produced with 
decreasing the altitude. Also, substantial differences in 
the overall kinetic temperature are observed when the 
axisymmetric case is compared to the 2-D case. For in- 
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Fig. 12 Kinetic translational temperature ratio (Tt /Too) 
profiles along the capsule surface for sections corresponding 
to s/R of (a) 0.4, (b) 0.8 and (c) 1.6. 



stance, the overall kinetic temperature core is located 
closer to the capsule nose for the axisymmetric case 
than that for the 2-D case, showing that the upstream 
disturbance due to the presence of the capsule becomes 
more pronounced for the 2-D case. 
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Fig. 13 Overall temperature ratio (To /Too) contours at the 
vicinity of the capsule nose for altitude of (a) 100 km, (b) 90 
km, and (c) 80 km. 



6.5 Computational and Experimental Comparisons 

Having completed the discussion on the primary prop- 
erties along the stagnation streamline and adjacent to 
the capsule surface, the attention now is turned to the 
comparison of the results. Nowadays, this is not a sim- 
ple task, given the small number of studies on SARA 
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Fig. 14 Comparison of wall pressure ratio (p w /Po) distribu- 
tion along the body surface in a (a) linear scale and in a (b) 
logarithm scale. 



capsule. As a result, the alternative is the comparison 
of the simulation data with those obtained experimen- 
tally with a similar geometry, i.e., a body defined by 
a spherical nose with a conical afterbody. In view of 
this difficult, Figs. 14(a-b) display the pressure acting 
on the body surface. In this set of figures, the wall pres- 
sure p w is normalized by the pressure at the stagnation 
point p , and the arc length s along the body surface is 
normalized by the nose radius R. In addition, Fig. 14(a) 
illustrates the pressure ratio p w /p distribution with a 
linear scale in the abscissa, while Fig. 14(b) depicts the 
same pressure ratio distribution in a logarithm scale in 
order to emphasize the pressure ratio behavior at the 
vicinity of the stagnation point. 

Based on this set of figures, Case A [35] refers to 
experimental data obtained in the hypervelocity tun- 
nel of the Arnold Engineering Development Center- von 
Karman Facility (AEDC-VKF) for a 9-degree half-angle 
spherically blunted cone. In addition, N2 was used as 
the working fluid, freestream Mach number near 19 
and freestream Reynolds number between 8,000 and 
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Table 4 Flow and geometric conditions for the spherically 
blunted cone models 



Case 


A [35] 


B [36] 


C [37] 


Moo 


~ 19 


~ 14.5 


5.8 


fleoo/in 


9 - 15 x 10 3 


5.7 - 12 x 10 3 


1.910 x 10 5 


R/ Rb 


0.3 


0.3 


0.8 


8 


9 


9 


10 



15,000/in. Case B [36] corresponds to experimental data 
obtained from Cornell Aeronautical Laboratory (CAL) 
shock tunnel for a 9-degree half-angle spherically blunted 
cone, with air as a working fluid, freestream Mach num- 
ber near 14.5 and freestream Reynolds number between 
5,700 and 12,000/in. Finally, Case C [37] represents a 
spherical nosed cone with semi- vertex angle of 10 de- 
grees tested in a hypersonic wind tunnel at a Mach 
number of 5.8, and freestream Reynolds number in the 
range of 0.97-2.38x 10 5 /in. Some characteristics of the 
models for Case A, B and C are tabulated in Table 4. 

According to Figs. 14(a-b), it is clearly seen that 
the wall pressure ratio for the present work agreed rea- 
sonably well with that obtained experimentally. It is 
important to remark that the shape of the pressure dis- 
tributions as described in nondimensional coordinates 
is independent of the spherical nose radius and of the 
Reynolds number over the range investigated. 



7 Concluding Remarks 

Computations of a rarefied hypersonic flow over the 
SARA capsule have been performed by using the Direct 
Simulation Monte Carlo method. The calculations pro- 
vided information concerning the nature of the flowficld 
structure about the primary properties at the vicinity 
of the nose and immediately adjacent to the afterbody 
surface by considering planar two-dimensional and ax- 
isymmetric geometries. 

Effects of rarefaction on the velocity, density, pres- 
sure, and temperature for a representative range of pa- 
rameters were investigated. The altitude varied from 
100 to 80 km, which corresponded to Knudsen num- 
bers Ktir from 0.4615 to 0.0115, Reynolds number Ren 
from 92 to 15249, and Mach number from 27 to 
29. Therefore, cases considered in this study covered 
the hypersonic flow in the transitional flow regime. 

It was found that changes on the altitude as well as 
on the capsule geometry disturbed the flowfield around 
the capsule, as expected. The domain of influence de- 
creased by decreasing the altitude. In addition, the do- 
main of influence for the axisymmetric case was smaller 



than that for 2-D case. Moreover, the extent of the flow- 
ficld disturbance along the stagnation streamline due to 
changes on the altitude was significantly different for 
each one of the primary flow properties. The analysis 
showed that the domain of influence for temperature 
is larger than that observed for pressure, density and 
velocity. 

The present document has described an initial inves- 
tigation of a high-altitude and low-density flow over the 
SARA capsule. Although this investigation has taken 
into account a representative range of altitudes, im- 
provements are still desirable to a realistic capsule de- 
sign. Since no database exist for such a design and since 
no appropriate flight data are available at such high en- 
try velocities and conditions, the use of computational 
methods is essential. Since many of these issues occur in 
the rarefied portion of the trajectory, accurate DSMC 
analysis are especially important. 
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Appendix 

This section focuses on the analysis of the influence of 
the cell size, the time step, and the number of molecules 
per computational cell on the surface properties in or- 
der to achieve grid independent solutions. Heat transfer, 
pressure and skin friction coefficients were selected in 
order to elucidate the requirements posed for the grid 
sensitivity study. As an illustrative example, the analy- 
sis presented in this section is limited to the capsule at 
an altitude of 90 km. The same procedure was employed 
for the other cases. 

A grid independence study was made with three dif- 
ferent structured meshes in each coordinate direction. 
The effect of altering the cell size in the ^-direction was 
investigated with grids of 75 (coarse), 150 (standard) 
and 225 (fine) cells, and 60 cells in the ^-direction. In 
addition, each grid was made up of non-uniform cell 
spacing in both directions. The effect of changing the 
cell size in the ^-direction is illustrated in Figs. 15(a- 
c) as it impacts the coefficients of heat transfer Ch 
(= 2q w /p 00 V^), pressure C p (= 2(p w - PooVPooO, 
and skin friction Cf (= 2T w /p 00 V^ ). The comparison 
shows that the calculated results are rather insensitive 
to the range of cell spacing considered. 
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Fig. 15 Effect of altering the cell size along the ^-direction 
on the coefficients of (a) heat transfer Ch, (b) pressure C p , 
and skin friction C/. 



Fig. 16 Effect of altering the cell size along the ^-direction 
on the coefficients of (a) heat transfer CV, (b) pressure C p , 
and skin friction C/. 



In analogous fashion, an examination was made in 
the 77-direction with grids of 30 (coarse) , 60 (standard) 
and 90 (fine) cells, and 150 cells in the ^-direction. From 
the total number of cells in the ^-direction, 60 cells are 
located along the spherical nose and 90 cells distributed 
along the conical afterbody surface. In addition, each 
grid was made up of non-uniform cell spacing in both di- 



rections. The sensitivity of the cell size variations in the 
77-direction on the heat transfer, pressure and skin fric- 
tion coefficients is displayed in Figs. 16(a-c). Results for 
the three grids arc basically the same, an indication that 
the standard grid was rather insensitive to the range of 
cell spacing considered. Therefore, the standard grid, 
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120 x 110 cells, for the 90 km case, is essentially an 
independent grid. 

One of the fundamental requirements in the DSMC 
method is that the time step At must be smaller than 
the local mean collision time. This requirement is re- 
lated to the important assumption employed in the 
DSMC method, i.e., the molecular movement can be de- 
coupled from the intermolecular interactions in a dilute 
gas if a sufficient small time step is used. In addition, in 
order to maintain a uniform distribution of simulated 
particles in the entire computational domain, a differ- 
ent time step At and scaling factor Fn can be obtained 
for each cell. Fn is the number of real particles rep- 
resented by one single simulated particle. As a result, 
the DSMC efficiency increases, and the computational 
effort is balanced within the simulated domain. It is 
worthwhile to highlight that although the time step At 
and scaling factor Fn vary among the cells, the ratio 
Fn I At must be the same in the entire domain. This re- 
quirement assures that the "mass flux" across the cell 
boundaries is conserved. 

With this perspective in mind, the following proce- 
dure is followed: (1) a computational grid is generated 
based on freestream conditions; (2) At and Fn values 
are defined for each cell according to the DSMC re- 
quirements and subject to the condition that Fn / At 
has the same value in every cell; (3) the parameters 
At and Fn are iteratively modified as the flow evolves 
within the simulated domain until each cell contains, on 
average, the desired number of simulated particles; (4) 
for the entire flowfield, all DSMC requirements are veri- 
fied, i.e., cell size smaller than the local mean free path, 
the time step smaller than the time related to the local 
collision frequency and a number of molecules around 
20-30 molecules. If within any cell these conditions are 
not satisfied the grid adaptation procedure, steps (1), 
(2) and (3), is restarted for a more appropriate spatial 
discretization. 

In doing so, At and Fn will be different for coarse, 
standard, and fine grids. By considering the standard 
grid for the 90 km case, a total of 13200 cells, the time 
step At changed from 4.5979 xl0~ 10 to 1.3953 xl0~ 5 , 
and the scaling factor Fn changed from 5.0601 xl0 +11 
to 1.5355xl0+ 16 . 

In a second stage of the grid independence investi- 
gation, a similar examination was made for the number 
of molecules. The standard grid for the 90 km case, 
120 x 110 cells, corresponds to, on average, a total 
of 189,100 molecules. Two new cases using the same 
grid were investigated. These two new cases correspond 
to 94,500 and 283,700 molecules in the entire compu- 
tational domain. The influence on the heat transfer, 
pressure and skin friction coefficients due to variations 
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Fig. 17 Effect of altering the number of molecules on the 
coefficients of (a) heat transfer Ch, (b) pressure C p , and skin 
friction C/. 



in the number of cells is depicted in Figs. 17(a-c). As 
the three cases presented approximately the same re- 
sults, hence the standard grid with a total of 189,100 
molecules is considered enough for the computation of 
the aerodynamic surface quantities. 

As part of the validation process, the axisymmet- 
ric version of this DSMC code was applied to a flat- 
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Fig. 18 Distribution of (a) normal velocity, (b) translational 
temperature, and (c) rotational temperature along the stag- 
nation streamline for the flat-ended cylinder. 



ended circular cylinder in a rarefied hypersonic flow. 
Results for velocity, translational temperature, and ro- 
tational temperature distributions along the stagnation 
streamline were presented and compared with those ob- 
tained from another established DSMC code [38] based 
in a code-to-code comparison. Since this comparison 
was published elsewhere [39], details will be kept to a 



minimum and the discussion restricted to the signifi- 
cant conclusions. 

In the computational solution, it was assumed that 
the flat-ended circular cylinder is immersed in a uniform 
stream flowing parallel to the cylinder itself. The cylin- 
der is modeled with a frontal-face radius R of 0.0185 
m, which corresponds to a 6Aoo, and a total length L 
of 30Aoo, where Aoo = 3.085 x 10~ 3 stands for the 
freestream mean free path. The freestream velocity Voo 
is assumed to be constant at 2694 m/s, which corre- 
sponds to a freestream Mach number M M of 10. The 
wall temperature T w is assumed constant at 570 K. 
This temperature is chosen in order to correspond to 
the temperature ratio T w /Too of 3.15 assumed by [38]. 

Normal velocity, translational and rotational tem- 
perature profiles along the stagnation streamline are 
illustrated in Figs. 18(a-c), respectively. In this set of 
plots, the normal velocity u is normalized by the veloc- 
ity Voo , the rotational Tr and translational Tt temper- 
atures are normalized by the freestream temperature 
Too. In addition, the distance x upstream the cylinder 
is normalized by the frontal-face radius R. Also, the 
solid curve represents the present DSMC simulations, 
and the full symbol represents the numerical data avail- 
able in the literature [38] . It is immediately evident from 
Fig. 18(a-c) that there is a close overall agreement be- 
tween both DSMC simulations at the vicinity of the 
stagnation region. 

References 

1. H. A. Carlson, J. of Spacecraft and Rockets, 36, 912 
(1999). 

2. M. R. Gilmore, R. Crowther, Proceedings of the 19th In- 
ternational Symposium on Rarefied Gas Dynamics, Eds. J. 
Harvey, G. Lord, (Oxford, University Press, 1995), 2, 1387. 

3. P. Gnoffo, Annual Review of Fluid Mechanics, 31, 459 
(1999). 

4. R. N. Gupta, J. N. Moss, J. M. Price, J. M., Proceedings 
of the 31st AIAA Thermophysics Conference, AIAA Paper 
1996-1859 (1995). 

5. M. S. Ivanov, G. N. Markelov, S. F. Gimelshein, L. V. 
Mishina, A. N. Krylov, and N. V. Grechko, N. V., J. of 
Spacecraft and Rockets, 35, 16 (1998). 

6. J. M. A. Longo, Aerospace Science and Techonology, 7, 
429 (2003). 

7. J. N. Moss, C. E. Glass, F. A. Greene, F. A., Proceed- 
ings of the 25th International Symposium on Rarefied Gas 
Dynamics, Saint Petersburg, Russia, (2006). 

8. R. Savino, M. S. Fumo, D. Paterna, and M. Serpico, M., 
Aerospace Science and Techonology, 9, 151 (2005). 

9. P. V. Vashchenkov, M. S. Ivanov, Proceedings of the In- 
ternational Conference on the Methods of Aerophysical Re- 
search, Russia, (2002). 

10. C. Weiland, J. Longo, A. Giilhan, and K. Decker, K., 
Aerospace Science and Techonology, 8, 101 (2004). 

11. R. G. Wilmoth, R. A. Mitcheltree, J. N. Moss, Proceed- 
ings of the 32nd AIAA Thermophysics Conference, AIAA 
Paper 1997-2510, (1997). 



18 



Santos, W.F.N. 



12. W. A. Wood, P. A. Gnoffo, and D. F. G. Rault, J. of 
Spacecraft and Rockets, 33, 643 (1996). 

13. F. Sharipov, Brazilian Journal of Physics, 33, 398 (2003). 

14. C. A. R. Pimentel, J. L. F. Azevedo, H. Korzenowski, M. 
B. H. Mantelli, Proceedings of the 43rd AIAA Aerospace 
Sciences Meeting and Exhibit, AIAA Paper 2005-0390, 
(2005). 

15. G. Tchuen, Y. Burtschell, and D. E. Zeitoun, Brazilian 
J. of Physics, 35, 148 (2005). 

16. P. G. P. Toro, M. A. S. Minucci, A. G. Ramos, J. B. 
Chanes Jr., A. L. Pereira, H. Korzenowski, H. T. Naga- 
matsu, L. N. Myrabo, Proceedings of the 39th AIAA 
Aerospace Sciences Meeting and Exhibit, AIAA Paper 
2001-0644, (2001). 

17. P. G. P. Toro, M. A. S. Minucci, J. B. Chanes Jr., A. G. 
Ramos, Proceedings of the 10th Brazilian Congress of Ther- 
mal Sciences and Engineering, Rio de Janeiro, RJ, Brazil, 
(2004). 

18. D. V. Kozak, F. Sharipov, Brazilian J. of Physics, 42, 
192 (2012). 

19. E. H. Hirschel, in Basic of Aerothermodynamics, Ed. 
by P. Zarchan, (Progress in Astronautics and Aeronautics, 
AIAA, 2005), 204, pp.17. 

20. G. A. Bird, in Molecular Gas Dynamics and the Direct 
Simulation of Gas Flows, (Oxford University Press, 1994). 

21. J. B. Pessoa Filho, Unpublished Communication, (2008). 

22. H. Machado, J. Boas, Proceedings of the 11th Brazilian 
Congress of Thermal Sciences and Engineering, Curitiba, 
PR, Brazil, (2006). 

23. C. Cercignani, in The Boltzmann Equation and Its Ap- 
plications, (Springer- Verlag, New York, 1988). 

24. G. A. Bird, G. A., 1981In Fisher, S. S., ed., Progress in 
Astronautics and Aeronautics: Rarefied gas Dynamics, Ed. 
by S. S. Fisher (AIAA, New York, 1981) 74, 239. 

25. G. A. Bird, in Rarefied Gas Dynamics: Theoretical and 
Computational Techniques, Ed. by E. P. Muntz, D. P. 
Weaver, D. H. Capbell, D. H., eds., (AIAA, New York, 
1989), 118, 374. 

26. C. Borgnakke, and P. S. Larsen, J. of Computational 
Physics, 18, 405 (1975). 

27. G. A. Bird, Physics of Fluids, 26, 3222 (1983). 

28. K. Guo, G. S. Liaw, in Proceedings of the 35th 
AIAA Thermophysics Conference, AIAA Paper 2001-2953, 
(2001). 

29. F. J. Alexander, A. L. Garcia, and B. J. Alder, Physics 
of Fluids, 10, 1540 (1998). 

30. F. J. Alexander, A. L. Garcia, and B. J. Alder, Physics 
of Fluids, 12, 731 (2000). 

31. A. L. Garcia, and W. Wagner, Physics of Fluids, 12, 2621 
(2000). 

32. N. G. Hadjiconstantinou, Physics of Fluids, 12, 2634 
(2000). 

33. R. F. Probstein, and N. H. Kemp, J. of the Aerospace 
Science, 27, 174 (1960). 

34. V. Cuda Jr., and J. N. Moss, J. of Thermophysics and 
Heat Transfer, 1, 97 (1987). 

35. B. J. Griffith, and C. H. Lewis, AIAA J., 2, 438 (1964). 

36. D. B. Wilkinson, S. A. Harrington, Rept. AF-1560-A-5, 
Cornell Aeronaut. Lab. 1962. 

37. R. M. Machell, in Aeronautical Engineer Thesis, CalTech, 
CA, (1956). 

38. D. I. Pullin, J. Davis, J. K. Harvey, in Rarefied Gas Dy- 
namics, ed. by J. L. Potter, (Astronautics and Aeronautics, 
AIAA, New York, 1977), 51, 379. 

39. W. F. N. Santos, in XII Encontro de Modelagem Com- 
putational, Rio de Janeiro, RJ, Brazil, (2009). 



